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1.  BACKGROUND  AND  INTRODUCTION 


Langmuir  (1948),  Bowen  (1950),  and  Ludlam  (1951)  early 
suggested  that  collision-coalescence  could  be  a  dominating 
mechanism  in  producing  rain  from  a  warmer  than  freezing  cloud. 

It  was  recognized  that  a  few  cloud  droplets,  after  having  grown 
to  a  size  larger  than  the  average  cloud  droplet,  would  fall 
faster,  and  also  sweep  out  cloud  droplets  more  efficiently. 
Findeisen  (1939)  had  earlier  considered  this  mechanism,  but 
doubted  its  efficiency.  He  and  Bergeron  (1935)  instead 
emphasized  mixed-phase  growth  (ice  and  water)  when  clouds 
extended  above  the  0°C  isotherm. 

By  calculating  values  for  the  collision  efficiency, 

Langmuir  showed  how  a  chain  reaction  could  initiate  precipitation 
after  seeding  with  large  water  drops.  Both  Bowen  and  Ludlam 
then  suggested  and  demonstrated  that  precipitation  could  result 
naturally  in  isolated  warm  convective  clouds.  Bowen,  however, 
concluded  that,  based  on  his  trajectory  calculations,  rain 
could  not  form  in  times  less  than  60  minutes.  He  realized 
that  this  conclusion  was  contrary  to  observation  in  that  most 
convective  clouds  go  through  a  growth-dissipation  cycle  in 
times  much  less  than  60  minutes.  Ludlam  concluded  that  collision- 
coalescence  is  probably  effective  only  in  maritime  clouds,  where 
there  is  more  opportunity  for  initial  formation  of  droplets 
greater  than  20  pm  radius. 

Because  rain  obviously  could  be  produced  in  warm  clouds  in 
times  less  than  an  hour,  other  methods  of  enhancing  the  coalescence 
process  were  suggested.  Cochet  (1951)  showed  how  droplet 
charging  could  provide  stimulation  to  the  collection  process, 
while  East  and  Marshall  (1954)  suggested  that  turbulence  could 
increase  collection  efficiencies. 

Recognizing  the  importance  of  collision  efficiency  values 
for  these  theoretical  studies,  Pearcey  and  Hill  (1956)  extended 
the  work  of  Langmuir  by  calculating  collision  efficiencies 
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for  droplets  of  comparable  size  up  to  200  ym  radius.  They 
neglected,  however,  interaction  of  each  droplet's  flow  field 
on  that  of  the  other.  Hocking  (1959)  concluded  that  such  an 
approximation  was  in  error  for  droplets  remaining  near  each 
other  for  a  period  of  time  (droplets  similar  in  size).  He 
removed  the  approximation  and  recalculated  efficiencies  for 
radii  up  to  30  ym.  His  results  differed  from  Pearcey  and 
Hill  in  that  he  found  a  limiting  upper  size  that  a  droplet 
can  collect,  whereas  the  former  found  collision  efficiency  to 
increase  as  the  ratio  of  the  two  radii  increased  to  one.  The 
most  significant  conclusion  of  Hocking,  however,  was  that  a 
droplet  of  radius  less  than  18  ym  could  not  be  a  collecting 
drop,  i.e.,  that  its  collision  efficiency  is  zero.  This 
conclusion  was  to  be  accepted  for  several  years. 

At  about  the  same  time  that  theoretical  studies  were 
initiated,  laboratory  investigations  to  determine  empirical 
values  for  collection  efficiencies  were  being  conducted.  Gunn 
and  Hitschfeld  (1951)  computed  growth  rates  for  1.59  mm  radius 
drops  falling  through  three  cloud  droplet  size  spectra.  They 
found  that  the  collection  efficiencies  calculated  from  these 
growth  rates  were  in  agreement  with  the  theoretical  work  of 
Langmuir.  They  concluded,  at  least  in  the  size  ranges  they 
considered,  that  coalescence  efficiency  is  equal  to  one.  Sartor 
(1954)  conducted  experiments  involving  small  droplets  of  similar 
size  (distilled  water  falling  in  mineral  oil).  He  found  that 
collision  efficiencies  for  sizes  from  5  to  15  ym  radius  were 
larger  than  those  calculated  by  Langmuir.  Telford  et  aj_. 

(1955)  also  attempted  to  determine  the  efficiencies  of  droplets 
nearly  equal  in  size.  They  found  that  such  droplets,  at  least 
in  the  size  range  75  ym  and  larger,  do  have  a  high  probability 
of  colliding  and  attributed  this  result  to  wake  effect  to  the 
rear  of  the  leading  drop. 

Schotl and  and  Kaplin  (1956)  conducted  experiments  using 
small  steel  spheres  falling  in  a  sugar  solution.  They  found 
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finite  values  of  collision  efficiency  for  sizes  of  5.5,  7.0, 

8.5,  and  10.5  pm  radius.  This  result  was  to  be  in  contradiction 
to  Hocking’s  theoretical  conclusions  two  years  later.  Hocking 
concluded,  however,  that  these  empirical  results  were  not 
representative  of  conditions  in  air.  Telford  and  Thorndike 
(1961)  gave  added  impetus  to  Hocking's  claim  when  their  results 
indicated  that  collisions  involving  droplets  of  about  35  ym 
diameter  were  very  rare. 

Further  refinements  to  theoretically  calculated  collision 
efficiencies  came  in  the  sixties.  Shafrir  and  Neiburger  (1963) 
extended  the  work  of  Hocking  by  calculating  efficiencies  for 
radii  up  to  136  ym.  Davis  and  Sartor  (1967)  provided  the  most 
significant  contribution  when,  after  making  further  refinements 
to  Hocking's  method  of  calculation,  they  found  that  there  was 
no  cutoff  below  18  ym  as  Hocking  had  found.  Collision  efficien¬ 
cies,  although  small,  were  finite.  This  was  a  significant 
finding. 

Although  theory  had  become  increasingly  refined,  laboratory 
work  gave  conflicting  results.  Collision  efficiencies  were 
not  always  verified.  As  several  of  the  early  researchers  had 
suggested,  however,  the  difficulty  might  lie  in  the  fact  that 
all  collisions  do  not  necessarily  result  in  coalescences,  i.e., 
coalescence  efficiency  is  not  equal  to  one.  Magono  and  Nakamura 
(1959)  photographed  collisions  between  a  large  stationary  water 
drop  and  smaller  colliding  droplets.  They  found  that  coalescence 
and  rupturing  of  the  larger  drop  is  a  function  of  the  smaller 
droplet's  size.  By  projecting  100-400  ym  radius  drops  on  a 
large  liquid  hemispherical  surface,  Schotland  (1960)  ascertained 
that  coalescence  was  related  to  two  dimensionless  parameters. 
These  parameters  are  a  function  of  drop  density,  medium  density, 
surface  tension,  normal  component  of  impact  velocity,  and  drop 
diameter.  By  impacting  water  droplets  with  radii  of  60  to 
200  ym  on  a  plane  water  surface,  Jayaratne  and  Mason  (1964) 
found  that  coalescence  is  a  function  of  the  angle  of  impact  and 
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that  larger  drops  coalesce  at  smaller  impact  angles. 

More  recently,  Whelpdale  and  List  (1971)  observed  collisions 
between  a  raindrop-size  drop  (500-1750  ym  radius)  and  a  droplet 
of  cloud  droplet  size  (35  ym  radius).  They  determined  that 
coalescence  occurs  for  85-95%  of  collisions  for  those  cases 
where  impact  velocity  was  equivalent  to  the  difference  in  their 
terminal  velocities.  As  others  had  also  found,  a  barrier  to 
coalescence  is  the  film  of  air  which  exists  between  the 
colliding  drops.  Based  on  their  observations,  they  found  that 
the  expression  E  =  r  ^  /  ( r ^  +  r  )  ,  where  E  is  the  coalescence 
efficiency,  r^  the  radius  of  the  large  drop,  and  rn  the  small 
droplet  radius,  is  valid  for  drop  radii  400  <  r..  <  2000  ym 
and  droplet  radii  20  <  r  <  100  ym. 

The  experiments  of  Whelpdale  and  List  concentrated  only 
on  rather  small  droplet/drop  (collected/collector)  ratios,  and 
then  only  for  rather  large  collector  drops.  Based  on  their  data, 
there  was  no  evidence  to  suggest  that  their  formula  was  really 
valid  for  increasing  size  ratios.  Levin  e_t  a_l_.  (  1973)  helped 
fill  this  gap.  Using  collector  drops  averaging  in  radius  from 
50  to  100  ym  and  the  collision  efficiency  data  of  Shafrir  and 
Neiburger  (1963),  they  computed  coalescence  efficiencies  for 
droplet/drop  ratios  of  0.18,  0.20,  0.22,  and  0.24.  Although 
their  values  fall  below  the  curve  suggested  by  Whelpdale  and 
List,  the  trend  is  the  same  --  decreasing  values  of  coalescence 
efficiency  with  increasing  size  ratios.  Feeling  that  this 
trend  is  quite  real,  Levin  ejt  aj_.  ,  suggested  that,  for  large 
size  ratios,  previous  results  in  which  collection  efficiencies 
agreed  with  computed  collision  efficiencies  probably  reflect 
enhanced  coalescence  through  droplet  charging  caused  by  the 
experimental  techniques  used. 

Since  the  continuous  growth  calculations  of  Bowen  (1950) 
and  Ludlam  (1951),  inclusion  of  the  stochasticity  of  the 
collection  process  in  numerical  computations  has  provided  an 
explanation  as  to  how  col  1 i s i on -coa 1 escence  can  lead  to  rain 
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within  20  to  30  minutes.  Telford  (1955)  was  the  first  to  make 
the  suggestion  that  a  discrete,  chain  reaction  type  process 
could  lead  to  rapid  growth  of  a  small  fraction  of  drops.  All 
earlier  growth  computations  had  assumed  continuous  growth 
whereby  all  droplets  of  a  given  initial  size  grew  equally 
with  time.  Extending  Telford's  work  which  considered  only 
two  initial  categories  of  drops,  Twomey  (1964)  studied  the 
evolution  of  a  spectrum  of  cloud  droplets.  In  general 
agreement  with  Telford,  he  found  that  production  of  a  few 
large  drops  proceeds  about  ten  times  faster  than  computations 
involving  continuous  growth.  Bartlett  (1966)  concentrated  on 
the  growth  region  from  20  to  40  ym  radius  and  determined  the 
effect  of  small  changes  in  the  collection  efficiencies  -- 
changes  which  were  to  simulate  electrical  or  turbulent  effects. 

He  found  that  any  external  alteration  which  produces  changes  of 
less  than  10%  can  be  ignored.  In  an  important  paper.  Berry 
(1967)  added  considerable  depth  to  the  study  of  droplet 
collection.  One  of  his  conclusions  was  that  continuous  growth 
is  important  initially,  before  a  second  spectrum  maximum 
of  water  mass  is  formed.  Thereafter,  stochastic  growth  becomes 
increasingly  important. 

New  computers  of  increased  speed  and  size  have  enabled 
researchers  to  utilize  collection  equations  in  cloud  models. 

Many  of  these  models  emphasized  the  dynamics  of  clouds  (e.g., 
Simpson  and  Wiggert,  1969;  and  Weinstein,  1970)  and  utilized 
parameterized  microphysics  as  developed  either  by  Kessler 
(1969)  or  Berry  (1968).  Others,  however,  did  include  explicit 
coalescence  growth:  Nelson  (1971),  for  the  purpose  of  developing 
rain  from  an  isolated  tropical  cumulus  with  a  time- independent 
updraft  profile;  Takeda  (1971),  in  a  two-dimensional  model 
which  concentrated  on  the  effect  of  an  environmental  vertical 
wind  profile;  and  Danielsen  et  al .  (1972),  for  the  purpose  of 
generating  hail.  Silverman  and  Glass  (1973)  showed  that  the 
results  of  parameterized  microphysics  in  cloud  models  are 
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significantly  different  from  those  involving  detailed  explicit 
microphysical  calculations. 

Implicit  in  all  of  the  above  cited  studies  is  the 
assumption  that  the  collection  efficiencies  are  identically 
equal  to  the  utilized  collision  efficiencies.  The  purpose  of 
this  paper  is  to  report  the  results  of  an  investigation  into 
the  effect  of  utilizing  the  empirical  data  of  Whelpdale  and 
List,  and  Levin  e_t  aj_.  for  the  coalescence  efficiencies.  The 
investigation  was  done  by  means  of  a  ti me- dependen t  eulerian 
model  of  stochastic  coalescence,  drop  breakup,  and  vertical 
mass  transfer.  By  utilizing  the  coalescence  efficiencies 
in  the  model  it  is  possible  to  determine  to  what  extent 
reduction  of  collection  efficiency  delays  the  onset  of 
precipitation  initiated  by  stochastic  coalescence,  an  onset 
that  can  be  demonstrated  to  occur  within  the  necessary  time 
frame  by  existing  theory,  but  a  theory  that  assumes  that  all 
collisions  result  in  coalescences.  If  the  delay  is  signif¬ 
icant,  efforts  should  be  directed  toward  other  positive  effects 
such  as  electrically  enhanced  coalescence.  At  the  least, 
it  should  suggest  continued  empirical  studies  of  coalescence 
efficiency. 
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2.  MODEL  DESCRIPTION  AND  NUMERICAL  PROCEDURE 


1.  BASIC  MODEL  DESCRIPTION 

The  basic  design  of  the  model  is  very  similar  to  that 
of  Nelson  (1971).  It  is  one -di mens i ona 1  with  ten  levels 
in  the  vertical.  The  purpose  of  the  model  is  not  to  simulate 
initial  cloud  development  nor  to  model  the  dynamic  aspects 
of  cloud  growth.  It  was  designed  for  studying  the  micro¬ 
physical  aspects  of  droplet  growth  in  a  tropical  warm  cumulus. 
For  this  purpose  a  droplet  spectrum  is  prescribed  initially 
at  each  of  the  model's  ten  levels  and  allowed  to  evolve 
by  stochastic  coalescence,  drop  breakup,  and  droplet  fallout. 
As  suggested  by  Nelson,  a  time- i ndependent  updraft  profile 
is  utilized  to  provide  a  way  of  accumulating  water  in  the 
upper  regions  of  the  cloud  and  to  provide  a  realistic 
restraint  for  hydrometeor  fallout. 

Forty-nine  categories  of  droplet  size  are  utilized  for 
spectrum  description  --  ranging  from  one  to  2820  pm  radius  in 
a  logarithmic  progression.  As  the  spectrum  evolves  in  time, 
droplets  enter  category  sizes  with  terminal  velocities 
exceeding  the  updraft,  and  exit  the  cloud  base  as  rain.  A 
special  array  is  utilized  to  store  this  fallout  in  order  to 
examine  its  characteri si ti cs . 

2.  THE  DIFFERENTIAL  EQUATION 


Let  N(r,z,t)  be  the  average  number  of  droplets  of  radius 
r,  at  a  height  z,  and  time  t  (definitions  of  symbols  are  given 
in  Table  1).  The  partial  differential  equation  which  then 
describes  the  processes  given  above  is 

- jpg - L=  [N(r,z,t)  ( w  ( z  )  -  u  (  r  ,  z  ) )  ] 

+  Q(r,z,t)  +  3 ( r  ,  z , t )  ,  (1) 
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Table  1.  Definition  of  Symbols 


b 

Constant  used  in  Golovin  collection  kernel  (1500  sec"1) 

B 

Partitioning  factor  for  "splitting"  scheme  (dimensionless) 

B 1 

Partitioning  factor  for  linear  "splitting"  of  droplet 
number  (dimensionless) 

B2 

Partitioning  factor  for  logarithmic  "splitting"  of 
droplet  mass  (dimensionless) 

E(i,n) 

Coalescence  Efficiency  (dimensionless) 

f (Inr) 

Number  of  droplets  as  a  function  of  Inr  (no.  m“3) 

1 » n 

Droplet  category  indicators 

M 

Total  number  of  droplet  categories 

No 

N(r,z,t) 

Initial  droplet  concent ra ti on  (no.  m~3) 

o 

Droplet  concentration  (no.  m  )  as  a  function  of 
r,z,t 

P(1,n) 

Collection  probability  of  an  ith  category  droplet 
collecting  an  nth  category  droplet 

Q 

Source/sink  term  symbolizing  collision-coalescence 

ri 

Radius  of  droplet  of  category  i  (cm) 

t 

Time  (sec) 

V  ( 1 ,  n ) 

Collection  kernel  as  a  function  of  the  ith  and  nth 
category  droplets  (cm3  sec~l) 

vo 

v  ( r ) 

Mean  volume  of  the  droplet  spectrum  at  t  =  0  (cm3) 

Volume  of  a  droplet  of  radius  r  (cm3) 

X 

Maximum  distance  from  the  collector  drop's  center 
(measured  perpendicular  to  the  line  of  fall) 
outside  of  which  no  collision  would  occur  (cm) 

X1’X2 

Mass  conversion  factors  used  in  coalescence 
equation  (dimensionless) 

Yc 

Linear  collision  efficiency  (dimensionless) 

z 

Height  ( m ) 

3 

Source/sink  term  symbolizing  drop  breakup 

Y 

Multiplicative  factor  for  logarithmic  droplet 
sequence 

ym 

Micron  (10“4  cm) 
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where  w(z)  is  the  updraft  speed,  u(r,z)  is  the  terminal 
velocity  of  a  droplet  of  radius  r  at  height  z,  Q  is  a 
source/sink  term  to  account  for  collision-coalescence, 
and  6  is  a  source/sink  term  for  drop  breakup. 

3.  DROPLET  CATEGORIES 

The  process  of  stochastic  coalescence,  operating  on  only 
a  small  initial  number  of  different  sized  drops,  by  definition, 
results  in  a  myriad  of  new  size  combinations  after  just  several 
time  steps.  It  is  physically  difficult  to  keep  track  of  all 
combinations  for  very  long.  As  a  result,  Eq.  (1)  is  solved 
at  each  of  the  ten  levels  for  a  sequence  of  fixed  radii.  For 
the  process  of  collision-coalescence  and  drop  breakup,  newly 
created  droplets  of  different  size  are  refit  back  into  the 
fixed  discrete  sizes.  The  fitting  procedure  will  be  described 
in  the  following  section. 

Because  of  the  range  of  droplet  sizes  (usually  one  ym  to 
several  thousand  ym  radius)  that  are  typically  covered  in  a 
microphysical  model,  care  must  be  taken  in  choosing  a  size 
sequence.  A  linear  sequence  of  discrete  sizes  is  not  practical. 
For  example,  even  a  five  ym  separation  results  in  200  categories 
for  a  1000  ym  maximum  size.  Ten  levels  thus  results  in  2000 
solutions  to  Eq .  (1)  in  one  time  step.  The  complexity  of  the 
coalescence  calculations  makes  such  a  sequence  impractical. 

A  more  workable  sequence  is  a  logarithmic  one  which  allows 
good  resolution  at  the  smallest  sizes  where  it  is  generally 
needed,  but  also  permits  inclusion  of  the  larger  precipitation- 
size  drops  if  they  should  be  needed.  The  discrete  sequence 
that  was  chosen  progresses  according  to  Eq.  (2). 

rj  =  Y  rj-i  (2) 

The  subscript  j  designates  droplet  category,  and  r,,  the  smallest 
radius,  is  1  ym. 


Eq.  (3)  describes  a  desirable  criterion  for  y . 

rj3  +  (yr-j)3  1  (yV)3  (3) 

Such  a  restriction  all ows  no  droplet  size  combinations  (resulting 
from  coalescence)  to  advance  past  one  category  greater  than  the 
largest  of  the  two  combining  droplets.  The  cube  is  significant 
because  liquid  water  is  proportional  to  radius  cubed. 

The  factor  y  can  be  determined  by  solving  Eq .  (3). 

Letting  r.  =  r,  =  1  pm,  we  have 

/  -,  \  3  3  6  ,  , 

( 1 )  +  Y  ±  y  •  ( 4 ) 

By  using  an  equality  for  a  solution  and  letting  y  =  A,  we 
have 

A2  -  A  -  1  =  0.  (5) 

Utilizing  the  quadratic  formula  produces  one  real  positive 
solution,  A  =  1.618.  The  cube  root  of  A  yields  y  =  1.174. 

This  value  of  y  produces  the  equality;  hence  a  correct  two- 
place  inequality  results  from  1.18. 

Consequently,  we  arrive  at  a  numerical  sequence  of  radii: 
1.00,  1.18,  1.39,  1.64,  etc.  Presently,  49  droplet  categories 
are  utilized  in  the  model.  This  number  permits  a  maximum  radius 
of  2820  pm,  the  size  at  which  drop  breakup  is  permitted  to  occur. 

4.  SOLUTION  OF  THE  DIFFERENTIAL  EQUATION 

By  recognizing  that  the  processes  described  in  Eq .  (1) 
are  independent  of  one  another  allows  each  process  to  be  solved 
separately  in  each  time  step  (for  each  droplet  category),  and 
then  summed  to  produce  the  total  change  in  N(r,z,t),  the  primary 
output  of  the  model.  Such  separation  is  desirable  because  of 
the  complexity  of  the  coalescence  computations.  An  explicit 
finite-difference  scheme  is  used  for  the  advective  term. 
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Because  Q  and  g  can  be  considered  as  constants  in  Eq .  (1), 
the  stability  criteria  for  advection  then  determines  total 
stabi 1 i ty . 

Time  derivatives  are  forward  differenced  while  the  "donor 
cell"  or  flux  method  of  advection  is  used  for  the  spatial 
derivative  (Thoman  and  Szewczyk,  1966).  This  form  of  the  space 
derivative  allows  for  mass  conservation  and  takes  the  form 

l.[N(w  -  u)]  _  Nz(wz  -  u z )  -  N z ” 1  ( w z ” 1  -  uz_1 ) 

3  Z  “  "  AZ  •  (6) 


5.  DROPLET  TERMINAL  VELOCITIES 


The  mathematical  approximations  of  Wobus  e_t  al  .  (1971) 
are  used  for  the  droplet  terminal  velocities.  These  approx¬ 
imations  allow  for  terminal  velocity  calculation  as  a  continuous 
function  of  droplet  radius;  they  are  in  excellent  agreement  with 
Stoke's  law  for  radii  smaller  than  50  ym,  and  with  the  data 
of  Gunn  and  Kinzer  (1949)  for  larger  drops  at  standard  air 
density  and  viscosity.  Because  the  terminal  velocities  in 
this  model  are  to  apply  at  heights  considerably  above  the 
ground,  the  effect  of  non-standard  atmospheric  conditions  is 
included.  Consequently,  the  velocity  of  a  droplet  changes 
slightly  from  cloud  base  to  cloud  top.  These  velocities  are 
calculated  once  at  the  start  of  the  program  and  stored  in  a 
two-dimensional  array  for  later  use.  The  collection  kernel 
(to  be  described  later)  is  also  dependent  upon  these  velocities. 

6.  COLLISION-COALESCENCE  OF  DROPLETS 

The  important  relevant  process  simulated  in  this  study  is 
collision-coalescence  --  the  term  Q  in  Eq.  (1).  For  a  given 
level  and  time,  the  change  in  droplet  number  due  to  stochastic 
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coalescence  is 


i  - 1  M 

Q ( r .  ) A t  =  - N ( r .  )  z  P(i,n)  -  z  N(r  )  P(n,i) 

1  1  n=l  n=i+l  n 

i  - 1 

+  N ( r . )  z  P  ( i  ,  n  )  X ,  ( i  ,  n  )  (1  -  B ( i  ,  n ) ) 

1  n=  1  1 

i  -2 

+  N ( r .  - )  z  P ( i - 1 , n )  X?(i  ,n)  B(i-l,n)  ,  (7) 

1-1  n= 1  L 


where  i  and  n  both  refer  to  category  numbers,  M  is  the  total 
number  of  categories,  P(i,n)  is  the  collection  probability, 

B  ( i  ,  n )  is  the  partitioning  factor,  and  Xj,  X2  are  mass  conver¬ 
sion  factors,  necessary  if  mass  is  partitioned. 

P(i,n)  is  the  collection  probability,  the  probability 
that  an  ith  droplet  will  collide  and  coalesce  with  an  nth 
droplet  in  time  At: 

P(i,n)  =  V(i,n)  E(i,n)  N ( r n ) A t  (8) 

where  V(i,n),  as  defined  by  Berry  (1967),  is  the  collection 
kernel : 

V(i,n)  =  Trri2Yc(ri  ,rn)2  (w(ri)  -  w(rn)).  (9) 

Yc  is  the  linear  collision  efficiency  used  by  Shafrir  and 
Neiburger  (1963).  As  defined  earlier,  u(r,z)  is  the  terminal 
velocity  of  a  droplet,  and  E(i,n)  is  the  coalescence 
efficiency  which  will  be  considered  separately. 

The  collection  kernel  describes  the  effective  rate  at 
which  volume  is  swept  out  by  the  collector  drop.  It  is  a 
function  of  the  area  through  which  the  collector  drop  falls 
(■nr.  ),  the  difference  between  its  own  terminal  velocity  and 
that  of  the  droplet  it  is  collecting  (u(ri5z)  -  u(r  ,z)),  and 
the  efficieny  with  which  droplets  are  swept  from  that  volume 
(Yc2).  This  collision  efficiency  is  defined  as: 


Y 


c 


2 


(10) 
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where  x  is  that  maximum  distance  from  the  collector  drop's 
center  (measured  perpendicular  to  the  direction  of  fall) 
outside  of  which  a  collision  would  not  occur  (Fletcher,  1962). 

The  partitioning  factor  B  determines  how  droplet  mass 
or  number  is  split  between  droplet  categories  as  it  is 
"advected"  into  larger  radii  categories  as  the  result  of 
collision-coalescence.  The  adoption  of  a  partitioning  scheme 
was  originally  stimulated  by  Bartlett  (1966)  and  is  similar 
to  that  of  Kovetz  and  Olund  (1969)  who  used 


B  ( i  ,  n  )  =  B 1  ( i  ,  n  ) 


(r.3  ♦  r„3)  -  r.3 


n 


ri  +  1 


3  -  r,3 


(id 


1+1  "  i 

for  partitioning  of  droplet  number.  Since,  in  this  case, 
droplet  number  is  being  red i s tr i bu ted  ,  the  conversion  factors 
x l  and  X2  are  unnecessary  and  are  set  equal  to  one.  The 
advantage  of  this  scheme  is  that  both  number  and  mass  are 
conserved.  Unfortunately,  it  leads  to  an  unrealistic  pseudo¬ 
diffusion  of  mass  such  that  the  solution  is  "fast."  Depending 
on  the  purpose,  such  diffusion  may  not  be  significant.  However, 
in  the  current  situation  where  the  time-development  of 
precipitation  is  the  critical  outcome,  the  diffusion  was 
considered  unacceptable.  Consequently,  a  more  accurate 
solution  was  sought. 


After  a  trial  and  error  examination  of  several  splitting 
schemes,  the  desired  accuracy  (to  be  demonstrated  in  the 
following  section)  was  found  with  a  logarithmic  partitioning 
of  droplet  mass,  instead  of  a  linear  partitioning  of  droplet 
number : 


B  ( i  ,  n  ) 


B2 ( i  ,n  ) 


£(i »n) 


-15- 


(12) 


where  C ( i » n )  =  exp[ln  r.3  +  Bl(i,n)  (In  r.  +  1 3  -  In  r.3)]  (13) 

and  B 1 ( i  , n )  is  that  defined  in  Eq.  (11).  By  recognizing 
that  r.j  +  j  =  yri  =  1.18  r^  ,  Eq.  (13)  reduces  to 


5  ( i  .  n ) 


r.3  ( y 3 ) B 1 ( i  , n ) 


r . 3(  1 . 643032  )  B 1  ^ 1  >n)  . 


(14) 


It  should  be  recognized  that  Eq.  (12)  represents  nothing  more 
than  the  logarithmic  position  (between  the  two  discrete  droolet 
sizes)  of  the  newly  created  droplet  mass,  rather  than  the 
linear  position  given  in  Eq.  (11). 

Because  mass  would  not  be  conserved  if  droplet  number 
were  partitioned  according  to  Eq.  (12),  mass  must  be  partitioned 
instead.  Consequently,  the  mass  conversion  factors  X1  and 
%2  are  necessary  for  droplet  number  recomputation: 


X 1  ( i  ,  n  ) 


r . 

l 


3 


X  2  ( i  ,  n  ) 


r 


i  -1 


3 


+ 
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Droplet  number  is  now  not  conserved,  but  as  will  be  shown  in 
the  following  section,  spectrum  mass  evolution  is  much  better. 

Equation  (7)  is  solved  for  category  numbers  i  =  3  to 
i  =  M - 1 .  Categories  1,  2,  and  M  are  special  cases  and  must  be 
considered  separately  since  they  fall  on  the  ends  of  the  distri¬ 
bution,  resulting  in  the  loss  of  several  of  the  terms  in  Eq.  (7). 
Eq.  (7)  allows  for  all  combinations  of  collisions  that  either 
increase  or  diminish  the  number  of  droplets  in  a  category.  The 
only  overriding  restriction  is  that  one  droplet  can  have  but 
one  collision  per  time  step.  The  first  term  in  Eq.  (7)  allows 
a  droplet  to  be  removed  because  of  a  collision  with  a  droplet 
smaller  than  itself.  The  second  term  allows  for  loss  due  to 
larger  collecting  drops.  The  third  term  allows  for  a 
droplet  to  re-enter  that  size  category  from  partitioning 
resulting  from  collisions  that  itself  incurs  (thus  producing 
a  droplet  size  that  must  be  partitioned  between  itself  and  the 
next  largest  category).  The  fourth  term  allows  for  re-entry 
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resulting  from  collisions  incurred  by  the  next  smallest  category 
(producing  a  droplet  of  a  size  between  itself  and  the  category 
in  question).  By  utilizing  the  criteria  imposed  by  Eq.  (3) 
there  is  no  danger  of  droplets  entering  a  category  from  sources 
other  than  those  described  above. 

7.  ACCURACY  OF  COALESCENCE  SCHEME 

In  order  to  test  the  accuracy  of  any  scheme  it  is  necessary 
to  compare  the  results  of  a  given  case  with  a  known  exact 
solution  of  that  same  case.  Golovin  (1963)  was  the  first  to 
derive  an  exact  analytic  solution  for  one  limited  case  --  one 
in  which  the  collection  kernel  was  proportional  to  the  droplet 
volumes: 


V  ( i  ,  n  )  =  b(j  it)  (  ri  3  +  rn3) 


(16) 


where  b  is  a  constant  (1500  sec-3).  Scott  (1968)  later 
provided  analytic  solutions  for  other  specific  kernels. 

Considering  both  linear  partitioning  of  droplet  number 
and  logarithmic  partitioning  of  droplet  mass,  the  above 
kernel  was  run  for  the  initial  droplet  number  spectrum 

O 

(normalized  to  1  gm  m"  )  used  by  Golovin  (here  on  a 
1 ogari thmi c  scale): 


f (lnr) 


3v  Nq 

— - —  exp ( -v/Vg) 


(17) 


where  v  is  the  volume  of  a  drop  of  radius  r,  Vg  is  the  mean 
volume  over  the  spectrum  at  t  =  0,  and  Mg  is  the  total  number 
of  drops  per  unit  volume  (also  at  t  =  0). 

In  Figure  1,  the  model  mass  spectrum  after  1800  sec  is 
compared  to  the  exact  analytic  solution  of  Golovin  for  both 
cases  of  redistribution.  It  is  clear  that  a  linear  parti¬ 
tioning  of  droplet  number  does  result  in  an  incorrect  widening 
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DISTRIBUTION  AT  t  =0  FOR 
ANALYTIC  AND  MODEL  SOLUTION 


—  ANALYTIC  SOLUTION  AT  t  =30  MIN- 
LINEAR  SPLITTING  OF  DROPLET  NUMBER 
liiaiiin  LOG.  SPLITTING  OF  DROPLET  MASS 


12  5  25  so 
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Figure  1.  Analytic  and  model  solutions  to  Golovin's  sum  of 
masses  collection  kernel'. 

of  the  distribution.  Logarithmic  partitioning  of  mass, 
however,  results  in  a  reasonably  accurate  final  solution. 

As  mentioned,  partitioning  of  mass  does  not  permit  number 
conservation.  However,  a  plot  of  the  number  spectra  reveals 
that  mass  partitioning  is  still  superior  to  linear  splitting. 

The  constant  kernel  of  Scott  (1968)  was  also  run  for 
the  initial  spectrum  given  in  Eq.  (17).  The  relative 
comparisons  (for  the  two  cases)  are  very  similar  to  the 
results  generated  by  the  Golovin  kernel,  with  logarithmic 
partitioning  of  mass  giving  a  better  solution.  On  the  basis 
of  these  results,  logarithmic  partitioning  of  mass  was 
chosen  for  the  "splitting"  procedure. 
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8.  COLLISION  AND  COALESCENCE  EFFICIENCIES 


Collision  efficiencies  for  radii  greater  than  28  pm 
are  taken  from  Berry's  (Berry,  1967)  numerical  fit  of  the 
Shafri r-Nei burger  data  (Shafrir  and  Neiburger,  1963). 
Efficiencies  for  radii  smaller  than  28  pm  are  based  on  the 
calculations  of  Davis  and  Sartor  (1967). 

The  coalescence  efficiencies  used  for  this  study  are 
based  on  the  empirical  work  of  Whelpdale  and  List  (1971), 
and  Levin  ejt  al_.  (1  973).  Whelpdale  and  List  observed  collisions 
between  two  drops,  one  of  approximately  35  pm  radius  and  the 
other  ranging  from  500  to  1750  pm  in  radius.  This  combination 
essentially  describes  hydrometeors  collecting  cloud  droplets. 
They  determined  that,  at  least  for  collector  drop  radii 
400  <  r..  <  2000  pm  and  collected  droplet  radii  20  <  r  n  <  100  pm, 


E  ( i  ,  n  )  = 


r . 


1.0 


( r .  +  r  ) 
v  i  n  ' 


(1  +  rn/ri ) 
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Levin  ejt  aj_.  concentrated  on  collector  droplets  from  50 
to  100  pm  in  radius,  with  larger  size  ratios.  Their  data, 
together  with  that  of  Whelpdale  and  List,  is  shown  in  Figure 
2.  Their  experimental  points  fall  below  the  Whel pdal e-Li st 
curve,  but  substantiate  the  downward  trend  with  increasing 
droplet/drop  ratio. 

Two  sets  of  coalescence  efficiencies  were  taken  from 
Figure  2  for  the  model  experiments.  The  first  (Case  II)  is  the 
Whel pdal e-L i s t  formulation,  assuming  that  it  applies  over  the 
entire  range  of  cloud  drop  sizes.  The  second  (Case  III)  is 
a  formulation  based  on  the  combined  data  of  Whelpdale  and  List, 
and  Levin  aj_.  (the  dashed  line  extension  in  Figure  2).  A 
third  order  least  squares  fit  produced  maximum  errors  of  less 
than  2%: 


E  ( i  ,  n  )  =  A ( 0 )  +  A ( 1 ) P  +  A(2)P2  +  A(3)P3, 

where  P  =  r  /r .  , 
n  l 
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Figure  2.  Experimentally  derived  values  of  coalescence 

efficiency  as  determined  by  Whelpdale  and  List  (1971) ,  and 
Levin  et_  al_.  (1973)  .  Solid  line  is  the  relation  suggested 
Levinet  aiS  and  List;  dashed  extension  is  that  suggested  by 


A ( 0 )  =  0.9986810817  , 

A(l)  =  -1.786534384, 

A(2)  =  -7.746074296, 
and  A ( 3 )  =  -1.149299369. 

It  is  recognized  that  these  data  are  incomplete  and  in 
need  of  further  verification.  The  qualitative  substantiation 
of  Whelpdale  and  List's  work  by  Levin  et  al.,  however,  is 
significant  for  the  following  reason.  It  indicates  that 
utilization  of  a  coalescence  efficiency  affects,  primarily, 
droplets  of  similar  size. 

Berry  and  Reinhardt  (1973),  in  their  continued  development 
of  improved  parameteri zati on  schemes,  recognized  that  the 


20 


growth  of  the  hydrometeor  phase  is  the  result  of  two  processes: 
accretion  (hydrometeors  collecting  cloud  droplets)  and  self¬ 
collection  (hydrometeors  collecting  hydrometeors).  While 
accretion  can  be  approximated  by  continuous  collection 
(average  growth),  Berry  and  Reinhardt  point  out  that  self¬ 
collection  must  be  described  by  stochastic  collection  and  can 
be  quite  important.  Consequently,  a  reduction  in  collection 
growth  as  the  result  of  incomplete  coalescences  primarily  affects 
sel f-col 1 ecti on .  It  was  the  specific  purpose  of  the  study 
reported  here  to  determine  how  the  reduction  of  that  component 
affects  the  initiation  of  precipitation  from  a  warm  cumulus. 

It  can  be  demonstrated  by  existing  theory  that  initiation  can 
occur  within  the  necessary  time  frame.  The  theory,  however, 
assumes  that  all  collisions  result  in  coalescences. 

9.  DROP  BREAKUP 

The  term  3  in  Eq .  (1)  accounts  for  drop  breakup.  As  used 
by  Nelson  (1971),  the  empirical  work  of  D'Albe  and  Hidayetulla 
(1955)  provide  the  data  on  which  drop  breakup  is  computed.  They 
conducted  experiments  with  large  (4. 5-6.0  mm  radius)  drops  in 
order  to  determine  the  number  and  size  of  the  breakup  fragments. 
Fragment  number  and  size  are  interpolated  for  their  data  in 
order  to  obtain  a  redistribution  of  liquid  water  for  category 
49  of  the  discrete  sequence. 

In  wind  tunnel  tests  Blanchard  (1950)  found  that  the 
size  at  which  drops  become  unstable  is  a  function  of  the 
degree  of  turbulence,  but  that  extreme  deformation  occurs 
for  drops  larger  than  2.5  mm  radius.  For  this  reason,  2820  ym 
(category  49)  was  chosen  as  the  largest  radius  in  the  discrete 
sequence.  A  drop  entering  this  category  is  broken  up  into 
approximately  16  fragments  and  redi stri buted  as  shown  by  the 
histogram  in  Figure  3. 
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DROP  RADIUS  (pm) 


Figure  3.  Redistribution  pattern  for  drop  breakup. 

10.  COMPUTATIONAL  PROCEDURE 

For  a  particular  experiment,  the  following  sequence  of 
major  events  is  followed  by  the  model: 

a)  Data  is  read  in  for  model  initiation:  cloud  droplet 
spectra,  Davis-Sartor  collision  efficiencies,  updraft  profile, 
as  well  as  other  bookkeeping  information. 

b)  Shafri r-Nei burger  collision  efficiencies,  droplet 
terminal  velocities,  and  coalescence  efficiencies  are  calculated 
and  stored  for  model  use.  Droplet  size  categories  are  set  up; 
the  input  cloud  droplet  spectra  are  fit  into  these  categories. 

c)  Collision-coalescence  operates  on  droplet  spectra. 

d)  Drops  entering  category  49  are  broken  into  fragments 
and  redistributed. 
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e)  Droplets  are  advected  between  model  levels. 

f)  Output  of  droplet  spectra,  liquid  water  contents, 
rain  fallout  at  cloud  base,  etc.  is  given  at  specified  time 
intervals  for  all  ten  levels  of  the  model. 

g)  Steps  (c)  through  (f)  are  repeated,  at  5  second  time 
steps,  until  the  final  time  limit  is  reached.  Experiments  to  be 
described  were  run  for  a  minimum  of  30  minutes  cloud  time. 
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3.  INITIAL  CLOUD  CONDITIONS 


For  comparison  and  consistency  purposes  the  same  initial 

droplet  spectrum  and  updraft  profile  as  used  by  Nelson  (1971) 

are  utilized.  The  data  of  Braham  e_t  aj_.  (  1957  )  for  tropical 

cumuli  with  no  echoes  are  used  for  the  initial  spectra  for 

all  levels  in  the  cloud.  A  liquid  water  content  of  0.8  gm  m”^ 

an  average  value  as  measured  by  Braham  e_t  a_l_.  ,  results  in  the 

water  mass  spectrum  given  in  Figure  4.  Initial  number  concen- 

-3 

tration  is  64  cm  .  Measured  average  droplet  concentrations 
_  3 

were  48  cm  for  the  non-echo  cumuli. 

As  in  Nelson,  a  Az  of  400  m  is  chosen.  In  this  model 
400  m  results  in  a  cloud  depth  of  3600  m.  As  mentioned 
previously,  an  additional  level  exists  below  cloud  base 
into  which  rainfall  is  accumulated. 

The  time-independent  updraft  profile  has  a  maximum 
speed  of  4  m  sec  ^  at  level  7  of  the  model,  decreasing  to 
zero  at  cloud  base  and  top.  The  profile  serves  two  purposes. 
It  leads  to  a  convergent  zone  above  the  profile  maximum  and 
thus  an  accumulation  of  water.  Although  the  mechanism  of 
accumulation  is  unrealistic  from  a  continuity  standpoint, 
liquid  water  increase  with  cloud  height  is  consistent  with 
increases  produced  as  the  result  of  the  condensation  process 
(assuming  that  the  adiabatic  value  is  not  exceeded)  and  with 
observations  (e.g.,  Warner,  1955).  Secondly,  the  updraft 
profile  provides  a  realistic  restraint  for  hydrometeor 
fall  out. 
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Figure  4.  Initial  cloud  droplet  mass  distribution-  used  in 
model  experiments  (LWC  =0.8  gm  m-^) . 
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4.  MODEL  RESULTS 


Three  experiments  were  conducted  to  determine  the 
relative  effect  of  using  the  coalescence  efficiencies.  Case 
I  is  a  control  case,  assuming  a  coalescence  efficiency  of  one. 

Case  II  assumes  coalescence  efficiencies  as  calculated  by  the 
Whel p d a  1  e-Li st  formula,  Eq.  (18),  and  Case  III  utilizes 
coalescence  efficiencies  taken  from  the  combined  data  of 
Whelpdale  and  List,  and  Levin  et  al . 

A  set  of  three-dimensional  diagrams  of  droplet  mass  spectral 
change  offers  a  convenient  method  of  observing  changes  in  the 
droplet  distribution  with  time  and  cloud  height.  Figure  5 
shows  spectral  change  for  Case  I,  the  control  case.  The  slant 
axis  (at  each  of  the  levels)  is  a  logarithmic  scale  of  droplet 
radius  while  each  of  the  intersections  of  a  droplet  spectrum 
with  the  horizontal  axis  represents  one  minute  of  time.  The 
vertical  axis  is  a  measure  of  liquid  water  content  (LWC)  under 
the  curve.  Consequently,  at  each  of  the  model's  ten  levels  a 
pictorial  description  of  the  droplet  spectral  evolution  is 
presented.  This  evolution  proceeds  as  the  result  of  droplet 
movement  between  model  levels  (updraft  plus  droplet  fall  velocity) 
and  movement  along  a  spectral  curve  ( co 11 i s i on -coal es cence 
pi  us  breakup ) . 

As  Nelson  (1971)  found,  temporary  accumulation  of  water 
at  the  top  of  the  cloud  leads  to  a  generation  zone  for  drops 
large  enough  to  fall  through  the  updraft  profile.  A  maximum 
LWC  of  2.8  gm  m“3  exists  at  level  10  at  about  7  minutes.  The 
significance  demonstrated  by  this  control  case  lies  in  the  fact 
that  local  temporary  accumulations  of  only  moderate  size  are 
enough  to  sufficiently  initiate  the  collision-coalescence 
mechanism,  leading  to  rain  exiting  cloud  base  in  about  15 
minutes  (see  Figure  5).  Such  moderately  large  LWC's  are 
consistent  with  observations.  Ackerman  (1959),  for  example, 
found  maximum  water  contents  of  over  3  gm  m  3  in  tropical  cumuli 
of  the  depth  simulated  here.  He  also  found  that  LWC's  were 
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Figure  5.  Droplet  mass  spectra  as  a  function  of  time  and 
height  for  Case  I  —  coalescence  efficiencies  equal  to  one. 
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successively  larger  with  increased  height. 

The  corresponding  spectral  diagrams  for  Cases  II  and  III, 
shown  in  Figures  6  and  7,  represent  increasing  limitations  on 
droplet  growth  because  of  incomplete  coalescences.  The  obvious 
change  to  be  noted  is  the  increasing  accumulation  of  water  at 
level  10,  especially  for  Case  III.  The  sharp  decrease  of 
coalescence  efficiency  (with  increasing  ratio),  as  shown  in 
Figure  2,  is  responsible  for  this  delayed  growth  in  Case  III. 

Figure  8  is  a  time  plot  of  LWC  for  level  10  for  Cases  II  and 
III,  together  with  the  control,  Case  I.  In  all  three  cases, 

LWC  rises  until  that  point  where  collision-coalescence  growth 
produces  drop  sizes  with  terminal  velocities  sufficient  to 
penetrate  the  updraft.  Simultaneously,  the  large  LWC,  which 
exists  primarily  in  relatively  small  droplets,  is  quickly 
depleted  by  growth  of  these  large  particles  (coalescence 
efficiencies  remain  close  to  one  for  small  size  ratios).  Case 
II  differs  only  slightly  from  the  control;  whereas  the  LWC  rises 
to  2.8  gm  m  for  Case  I,  a  maximum  of  3.6  gm  m  is  reached 
for  Case  II.  Spectral  evolution  at  cloud  top  (for  Cases  I  and 
II)  is  a  continuous  process,  the  second  water  maximum  forming  as 
an  evolution  from  the  primary  maximum.  This  evolution  is 
delayed  two  to  three  minutes  as  the  result  of  the  use  of  the 
Whel pdal e-Li st  formula. 

Case  III,  however,  because  of  its  severe  limitation  on 
growth  between  drops  of  similar  size,  presents  a  very  different 
situation.  The  LWC  at  level  10  rises  to  more  than  double  that 
of  the  control.  Development  of  large  hydrometeors ,  when  it  does 
happen,  does  not  progress  as  a  continuous  evolution  from  the  bulk 
of  liquid  water.  At  about  20  minutes,  large  hydrometeors  appear 
di scontinuously  when  viewed  on  an  LWC  plot  (as  given  in  Figure  7). 

Rainfall  rate  and  cumulative  rainfall  are  plotted  in 
Figures  9  and  10.  One  might  initially  suspect  that,  although 
precipitation  is  delyaed  for  Cases  II  and  III,  rainfall  rate 
and  total  rainfall  would  be  as  large.  Such  is  not  the  case. 
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Figures  6.  Droplet  mass  spectra  as  a  function  of  time  and  height 

whJ-iSqf?,,1?  •  I  £oaleifcence  efficiencies  as  defined  by  the 
whelpdale-List  formula. 
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Figure  7.  Droplet  mass  spectra  as  a  function  of  time  and 
height  for  Case  III  —  coalescence  efficiencies  as  defined 
by  the  Whelpdale-List ,  Levin  et_  al_.  combined  data. 
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Figure  8.  Liquid  water  content  change  at  cloud  top  for 
Cases  I,  II,  and  III. 


Figure  9.  Rainfall  rates  for  Cases  I,  II,  and  III. 
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Figure  10 .  Cumulative  rainfall  for  Cases  I,  II,  and  III. 

Maximal  rates  of  approximately  18,  11,  and  9  mm  hr"1  occur  at 
20,  25,  and  32  minutes  for  Cases  I,  II,  and  III,  respecti vely . 
Similarly,  cumulative  rainfall  is  less  for  periods  up  to  40 
minutes. 

Raindrop  size  spectra  for  the  three  cases  were  compared 
to  each  other  and  against  the  Marshal  1 -Palmer  spectra  (Marshall 
and  Palmer,  1948).  For  rainfall  rates  common  to  the  three  cases 
it  was  found  that  the  spectra  differed  only  slightly  from  each 
other.  In  general,  however,  the  spectra  for  Cases  II  and  III 
were  shifted  toward  smaller  sizes.  For  example,  Figure  11  shows 
comparative  spectra  at  the  rainfall  rate  which  was  maximum  for 

Case  III,  9.5  mm  hr"  .  The  diagonal  line  is  the  Marshall -Palmer 
spectrum. 
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Figure  11.  Raindrop  spectra  just  above 

rainfall  rate  of  9.5  mm  hr“l  for  Cases  I,  II,  and  III” 
Diagonal  line  is  the  Marshall-Palmer  spectrum  for  that 

^  "t  0  • 
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5.  SUMMARY  AND  CONCLUSIONS 


The  stochastic  theory  of  collision-coalescence  has 
enabled  cloud  physics  to  explain  production  of  rain  from  a 
warmer  than  freezing  cloud  within  20  minutes.  Calculations 
which  can  demonstrate  this  occurence,  however,  are  based  on 
the  assumption  that  all  collisions  result  in  coalescences. 

Recent  empirical  studies  have  determined  that  coalescence 
efficiency  decreases  with  increasing  size  ratio  of  the  two  drops. 
A  time-dependent  one-dimensional  model  which  incorporates  the 
effects  of  stochastic  coalescence,  drop  breakup,  and  vertical 
transfer  by  means  of  a  constant  updraft  profile  was  used  here 
to  study  the  effect  of  using  these  efficiencies.  Utilizing 
two  separate  sets  of  data  which  outline  the  suggested  limits 
of  empirical  findings,  comparisons  were  made  to  a  control  case 
in  which  rain  could  be  generated  from  a  tropical  cumulus  in 
about  15  minutes. 

The  results  were  found  to  be  quite  sensitive  to  coalescence 
efficiencies  for  collisions  between  droplets  very  close  in  size, 
the  region  in  which  the  two  sets  of  data  differed.  In  general, 
however,  utilization  of  these  efficiencies  leads  to  delay  of 
precipitation  and  to  decreases  both  in  rainfall  intensity  and 
cumulative  rainfall.  For  equivalent  rates  of  rainfall,  however, 
raindrop  size  spectra  were  found  to  differ  only  slightly.  The 
most  significant  finding  is  the  large  LWC  which  is  necessary 
to  sufficiently  stimulate  the  collision-coalescence  mechanism 
for  those  cases  where  coalescences  are  restricted.  Using  the 
more  restrictive  set  of  coalescence  efficiency  data  resulted  in 
maximum  LWC 1  s  which  are  unreasonably  large  and  uncommon  for 
typical  tropical  cumuli.  If  coalescence  efficiencies  are  indeed 
as  limited  as  this  one  set  of  empirical  data  indicate,  then 
other  compensating  effects  (e.g.,  electrical)  must  be  found. 

It  is  suggested  that  more  laboratory  work  be  conducted  to  verify 
and  extend  the  empirical  results  thus  far  determined. 
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